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TECHNICAL NOTE 3964 


THE LINEARIZED SUBSONIC FLOW ABOUT SYMMETRICAL 
NONLIFTING WING-BODY COMBINATIONS 
By John B. McDevitt . 


SUMMARY 


Methods are presented for determ i ning the linearized subsonic flow 
about symmetrical, nonlifting wing-body combinations. The wing is repre- 
sented by a suitable planar distribution of elementary sources in accord- 
ance with the usual assumptions of linearized, thin-airfoil theory. The 
required boundary condition for tangential flow at the body surface is 
met by distributing along the body axis suitable distributions of three- 
dimensional sources and multlpoles. As part of the present analysis, a 
theory is presented for determining the flow about circular or noncircular 
bodies by the use of axially distributed sources and multlpoles. The flow 
over thin, nonlifting airfoils having symmetrical profiles is also studied 
in considerable detail, and it is shown that the concept of the oblique 
source filament can be used in obtaining numerical results for wings having 
tapered as well as untapered plan forms. 

A comparison of theory and experiment is made for two 45° sveptback 
wings in combination with basic Sears -Haack bodies of revolution and in 
combination with the basic bodies indented according to Whitcomb's tran- 
sonic area rule. The effect of indenting the body on the pressure dis- 
tribution near the wing-body juncture is of particular interest for swept- 
wing and body combinations. The effect of the area -rule indentation was 
to adjust the pressure field near the body so that the wing surface Isobars 
tended to follow the local sweep lines of the wing. 


INTRODUCTION 


Linearized theories for predicting the subsonic pressure distributions 
on thin wings or slender bodies have been developed in considerable detail 
but comparatively little attention has been directed toward the problem 
of predicting the subsonic pressure distributions on wings and bodies in 
combination. The existing theories for determining the subsonic flow about 
wing-body combinations are based on the assumption of extremely slender 
wings and bodies in accordance with the basic concepts of slender-body 
theory as initiated by Max Munk and R. T. Jones (refs. 1 and 2) and 
extended by numerous authors (see, in particular, refs. 3 to 6). 
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In the present paper the fundamental solutions of linearized, three- 
dimensional-flow theory are used in the development of a practical method 
for calculating the subsonic flow about nonlifting wing-body combinations. 
(The use of linearized theory requires, of course, that the wing be thin 
and the body slender but does not require the assumption of extreme slen- 
derness involved in slender-body theory. ) In the present method partic- 
ular emphasis is placed on the practical application of the theory. No 
attempt is made to estimate the accuracy of calculated results using this 
method but a comparison of theory and experiment is provided. 

In the analysis which follows,' general relationships are first derived 
for determining the strengths of the wing planar sources and the singular- 
ities at the body center line so as to satisfy the required boundary con- 
dition for tangential flow at the wing-body surface. Next, the body alone 
and the wing alone are considered, and finally, the wing and body in 
combination are considered. 


PRIMARY SYMBOLS 


A wing aspect ratio 

Anj^n strengths of axial raultipoles according to linearized theory 

b wing span 




pressure coefficient, local pressure minua free-stream static 
pressure divided by free-stream dynamic pressure 


c local wing chord 

c. reference chord near wing-body junction (chord through the 

J point of intersection of body and wing leading edge) 


°R 


chord at wing center line 


c^ chord at wing tip 


E 


f 




source strength 

body fineness ratio (body length divided by maximum body 
diameter) 

multipole strength parameter (see eq. (34)) 


A 









NACA TU 3964 


3 


l 

L 

Moo 

N 

r 

R 

Ra 

R 

AR 


plan-form convergence factor. — 

b/2 

"body length, (distance from nose to theoretical point of 
closure) 

lateral distance to extended wing vertex for p = 1 (see 
sketch (l)) 

free -stream Mach number 

see sketches (e), (f), and (m) 

see sketch (c) 

polar coordinate in yz plane 
body radius 

body radius at wing leading edge 

radius of equivalent circular body having equal cross- 
sectional area 

incremental distortion of body shape from the circular shape 
having equal cross-sectional area 

cross-sectional area of body 


t 

c 


wing thickness to chord ratio 


free-stream velocity 

u,v,w 

v (x)' v (e)' v (r) 

n, v,v 1 perturbation velocities normalized by division by the 

V (x)' V (y)’ V (z) J free-stream velocity 


j- total velocities 


V (r) 

V (t)' V (n) 


radial component of perturbation velocity in yz plane 

tangential and normal perturbation velocities for an oblique 
source filament (see sketch (f)) 
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azimuthal component of perturbation velocity in yz plane 

radial-perturbation velocity at body surface induced by wing 
sources 


A v/ 


v (r) v ' n r ) w J 


X. 


■1 

x,y,z 

Z 


'v i- - - "-av 
x - x 

Cartesian coordinates 
wing profile ordinate 


Z* 


> 


dZ 

dx 




5 


J 


I 




P 


7 

0 

A 


9 

A 

Ae 


airfoil slope in streamwise direction 
modified airfoil slope (see eq. (59)) 

Cartesian coordinates normalized by division by a reference 
chord 

Cartesian coordinate measured from the leading edge of the 
reference chord Cj (see sketch (o)) 

variable of integration (see sketches (g) and (l)) 

|-Pr) tan Ae 

JW-A 

s in 2 A+0 2 co s 2 A = *7 l-l^cos 2 A 

ratio of specific heats (for air y =1.4) 
polar angle in yz plane 

Crjn 

wing plan-form taper ratio , •— 

°R 

perturbation potential 

angle of sweep, positive when swept back 

effective angle of sweep according to the GOthert 
compressibility rule 
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( ) variable of integration 

( ) ' first derivative with, respect to the free-stream direction 

( ) " second derivative with respect to the free-stream direction 

multipole influence functions (see eqs. (28) to (30)) 

Subscripts 


av average 

b body 

e equivalent configuration of the Gothert compressibility rule 

j condition along wing-body junction 

LE wing leading edge 

TE wing trailing edge 

n order of multipole 

w wing 

y ; Zj0 partial derivatives when not enclosed by parentheses 

GENERAL THEORY 
Statement of the Problem 

In this paper methods are presented 
for determining the subsonic flow about 
nonlifting wing-body combinations. 

These methods are based on the usual 
assumptions of linearized flow theory. 

The flow over the wing is to be deter- 
mined by representing the nonlifting 
wing by a suitable distribution of 
planar sources. The wing is assumed 
to extend through the body as indicated 
in sketch (a). 



£ 

Sketch (a) 



6 


MCA TN 3964 


It is further assumed that the wing-body combinations possess both 
lateral and vertical symmetry. The required boundary condition for tan- 
gential flow at the body surface is to be satisfied by placing along the 
body axis (which is coincident with the wind direction) suitable axial 
distributions of sources and multipoles. 


Boundary Conditions for Tangential Flow 
at the Wing or Body Surface 


If F(x,y,z) = 0 is the general equation determining a body or wing 
shape in Cartesian coordinates , the boundary-condition for tangential 
flow at the surface is 


U — + V 

Sx 


— + w — 

by 5z 


= 0 


( 1 ) 


In cylindrical coordinates the boundary condition can be written as 

V, v— + V, — = 0 

(x) 5x (r)^ R B ^ 

where F(x,R,0) = 0 determines the body shape. 


( 2 ) 


Wing theory .- If the ordinates for a thin airfoil are represented 
by the function Z = Z(x,y), then 

F(x,y,Z) = Z-Z(x,y) = 0 

and the boundary condition at the surface is 


U — + - w = 0 

bx by 


( 3 ) 


For linearized, thin -airfoil theory, the second term can be neglected, 
U = U M (l+u) is replaced by Uco, 


dZ 

dx 


* 


w 


z — > 0 


( 4 ) 
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by Uoo and the boundary condition in terms 
of perturbation "velocities at the body surface 
is 

Sr v (e) Sr 

Sx B Se M 


Sketch (b) 


(6) 


Compressibility Corrections 


The rules of correspondence between linearized compressible and 
incompressible flow in two dimensions were given originally by Prandtl. 
The extension of these rules to include flows in three dimensions was 
made by Gothert (ref. 7). The use of the Gdthert rule is as follows: 
First, the lateral and vertical dimensions of the configuration are 
reduced by the factor 3 so as to obtain an "equivalent" configuration. 
The equivalent sweepback A e of the wing reference line (or oblique 
source filament) is obtained from 


tan Ag = 


tan A 


Thus, the quantities sin A and cos A are replaced by 


sin Ag = 


sin A 


cos Ag = 


*7 sin 2 A+3 2 cos 2 A 
3 cos A 

N/sin a A+3 2 co s 2 A 


sin A 

Pa 


3 cos A 

Pa 


(7) 


(8) 
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Next, the incompressible flow about the equivalent configuration is 
determined. Finally, the compressible flow is obtained from the calcu- 
lated incompressible flow by increasing the velocity potential by the 
factor l/P 2 , the axial -perturbation velocity is increased by the factor 
l/(3 2 , and the lateral and vertical velocities are increased by the factor 
l/p (the pressure coefficient is increased by the factor l/j3 2 ). 


Pressure Coefficient 


The exact expression for the pressure coefficient (isentropic flow) 
is 



(9) 


In order to obtain an approximation consistent with linearized theory this 
equation may be expanded to give (see, e.g. , ref. 8) 


Cp = -2u-(v 2 +w a )+p 2 u 2 +2^u(v e +w 2 ) + . . . (10) 


For the linearized flow past planar systems,- the first term is sufficient. 
For slender bodies it is necessary to consider additional terms and the 
appropriate expression is (in cylindrical coordinates) 


Cp = ~2u 



(H) 


Nonlifting Thin-Airfoil Theory 

The basic linearized partial differential equation for subsonic flow 
is the well-known Prandtl-Glauert equation 

(l-^^xx^yy^zz = 0 ( 12 ) 


where Moo is the free-stream Mach number and q>(x,y, z) is the perturbation 
potential. 
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Perturbation potential and, velocities . - The perturbation potential 
subject to the boundary condition of equation ( 4 ) can be written in the 
form (see ref. 9 for the fundamentals of source-sink theory) 


tp(x,y,z) 


-1 ff PE (5c, py) as dy 


( 13 ) 


where the region of integration £ extends over the entire plan form of 
the wing. This solution represents the potential due to a planar (z => 0) 
distribution of sources having strength E per unit area. 

The streamwise, lateral, and vertical perturbation velocities induced 
by the planar distribution of sources are: 


u ( ) = = 1 fT g(x-x)E(x,py)dx dy 

dx ft- ^ r ^2 . .2 jV 


2 „2r/ 2 


“ s -|(x-x) +3 [ (y-y) +z ] 


(HO 


v(x,y) = 


ckp 

Sy f 


g 2 (y-y)E(x,Py)dx dy 


ff 

S |(x-x) 2 +p 2 [ (y-7) 2 +z 2 ]j- 7 


( 15 ) 


/ n _ ^9 1 (T p g zE(x,3y)dx dy 

^ ' Sz " J . 2 r . _. 2 2l 1 3/2 

£ i(x-x) +3 [ (y-y) +z ] i 


(16) 


Determination of the source strength distribution . - An approximation 
for the source strength E can be obtained by finding a boundary value 
for w(x,y) as z -> 0 and using that result in conjunction with the 
boundary condition of equation ( 4 ). 

The boundary value for w as z — > 0+ is to be obtained from 


v = lim 

Z= ° + z^o+ 


I 


P 2 zE(x,3y)dx dy 


|(x-x) 2 


(y-y^+z* 3 


^3/2 

] } 
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09 t J8(y+dy) f 0y 


The above Integrand Is zero 
everywhere for z = 0 except 
at the point x = x, y = y 
where a singularity of the 
kind 0/0 occurs; conse- 
quently E(x,y) may he 
replaced by E(x,y) and 
removed from under the inte- 
-gral signs. In order to 
evaluate this equation it is 
convenient to make the fol- 
lowing changes of variables 
(see sketch (e)): 




sin -1 


P-7(y-y) 2 +z 2 

Ni 


^2 


sin" 1 


Pz 

N 2 


Sketch (c) 


where 


Nx 


J (x-x) 2 +p 2 [ (y-y) 2 +z^] 


p7(y-y) g + z2 

sin A-l 


n 2 = P N /(y-y0 2 +z 2 = 


ftz 

sin Ag 


Finally there is obtained. 


. f\ [* si » ^ 

4irPU«. J J 

E(x,py) 

2PU00 
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This result is subject to the "boundary condition (4) and. hence the source 
strength is defined as 


E(x,(3y) = aZg(x,Py) 

U ra dx 


(IT) 


■where dZy/dx is the airfoil slope measured streamwise along the upper 
surface of the wing. 

So lire e filaments . - In many cases, where the wing plan-form edges are 
straight-line elements, the use of source filaments in the mathematical 
representation of the flow is advantageous. In this way the determination 
of the perturbation potential is reduced from a two-dimensional integra- 
tion over the entire plan form to a single integration performed in the 
streamwise direction. Details concerning the use of source filaments 
will he presented later in this report. 


The General Solution in Cylindrical Coordinates for a 
Line of Axially Distributed Singularities 


The linearized potential equation in cylindrical coordinates is 

(l-lO'P xx +<Prr + J %! = 0 < l8 > 

The linearized spatial flow about bodies or wing-body combinations requires 
a consideration of suitable distributions of sources and multipoles along 
the body axis with the axis of the nrultipoles aligned with the wind axis. 
The general solution of equation (l8) consists of an infinite set of basic 
singularities referred to as nrultipoles (ref. 9, p. 527 )* 

Perturbation potential and vel ocities . - The perturbation velocity 
potential of a line of axially distributed nrultipoles of order n in 
compressible flow is 


<p( x ,r,0) = 




'I 


(0r) 


n 


cos 


no/ 


A n (x)dx 


(Pr) sin n© 


/ 


[(x-x) 2 +(pr) 2 H l/2 ) 
B n (x)dx 


iu-*) 2 + (M 2 r- (l/2) 

where A Q (x) and B^x) represent multipole strengths. 


( 19 ) 
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The perturbation velocity components can be obtained from equa- 
tion ( 19 ) by partial differentiation in the axial, radial, and azimuthal 
directions. If it is assumed that A 11 (0) = A n (l) = 0 and 1^(0) = B^Z) = 0, 
the velocity components are: 


u(x,r,0) = — ~y (2n+l)(pr) n f 
4*f3 2 Ls J 


^ (x-x)[A tl (x)cos n6+B n (x)sin n8] 
dx 

t(x-x) 2 + Or) 2 ] n+ < 3/2 > 


( 20 ) 


V( r )(^ r ^) 



[A n (x)cos n0+B n (^)sin n0] 
[(x-x) 2 +(pr) 2 ] n+(l/2) 


dX + 


1 

4rt0 


I 


(2n+l)(pr) n+1 



[A n (x)cos nO+BnfxJsin n0] 
[(x-*) 2 + Orff + < 3,2 > 


dx 


( 21 ) 


V( Q )(x,r,0) 



[A n (x)sin n0-B n (x)cos n0] 
[(x-S) 2 +(|3 3 -) 2 ] n+ ( l/2 ) 


dx 


( 22 ) 


If the change of variable 


T x-x 

5 = 17 

is made, the above equations can be written in the following forms 

Z-x 


(23) 


cp(x, 


,r,0) = | 1 . n / [A n (x+|3r5)cos n0 

4rtf3 2 U (pr) n 


-x 

pr 


B^x+prSjsin n0]f n (s)dS 


(24) 
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I-x 


u(x,r, 0 ) = 2n ’~r f Pr (A^cos ne+BaSln neJ^Cs)^ (25) 

W 2 A. (pr) + 

Sr 


I-x 


v ( r )( x , r , 0 ) = 


JLV — / (A n cos ne+B^sin n0) [ (2n+l)h n (s)-nf a (8) ]d£ 

^V(P r ) n+lv k 

Pr 


(26) 


v (0)( x,r,fl ) 


-JlY. 


I-x 


n 


^P 40 r) n+:L ^ 

Pr 


^ (A n sin nS-BaCos n0)%(s)d8 


n 


(27) 


■where 


f “ (5) ' (5W + < l/2 > 


(28) 


8^(5) = 


(& 2 +i) n+ ( 3/a) 


(29) 


h n (s) = 7— — r 

(B S +l) n+ 3/2 


(30) 


The variations of the influence functions f^Cs ) , ^(8), and 1^(8) 
with 8 are shown in figure 1 for values of n equal to 0, 1, and 2. 

An alternative expression for the axial perturbation velocity 
involving the influence function f n (s) is obtained from equation (25) 
by integrating once by parts 

___ I-x 

u(x,r,0) = — - I Pr (A n 'cos ne+Bn’sin n0)f n (8)d& 

^ ^ (.fr) J -x 


( 31 ) 
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where 


An' 


(Br 


dA n (x+prS') 

dX ™~ 


V 


dB n (x+3r$’) 

dX 


In the preceding results, the solution for an axial distribution of 
sources results when n is set equal to zero. For n = 1 the "dipole" 
solution appears; for n = 2, the "quadripole" solution; etc. For small 
values of 0r the attenuation of the induced perturbation velocities 
with increasing radial distance pr varies approximately as l/(Pr) n+1 . 

Limiting values for v ( r ) v (g) as P r 0.- ® le ( 3i rec '4 1136 

the preceding equations in boundary-value problems is difficult since the 
unknown multipole strength functions must be related to the boundary con- 
dition (6). An approximation for the strength functions in terms of 
boundary values can be obtained by first investigating the limiting values 
for v ( r ) and. v (g) as 0r-> 0. The following change of variable is made 

(see sketch (d)). 



A = sin" 


Pr 

7(x-x) a + p 2 r e 


As pr -> 0 it is permissible 
to replace A^x) and %(?) by 
A n (x) and Bn(x). Consequently 
there are obtained the follow- 
ing limiting values for 


r (r) 


and v 


( 0 ) 


as Pr -> 0: 


Sketch (d) 


n 


V) =I T w n ‘I 


K n 


4rtP(Pr) 


n+i 


[A n (x)cos n0+B a (x)sin n0] 


(32) 
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= ) — [A r ,(x)sin n0-B n (x)cos n0] 

n 4 M*(Pr ) n+1 


(33) 


■where 


Kq = K x - 2 1 


(g)(4) ■ . ■ (2n) 
(1)(3) ■ • . (2n-l) ^ 


(34) 


Equations (32) and (33) will "be used with appropriate boundary 
conditions in subsequent sections of this report in order to determine 
multipole strength distributions for slender bodies or wing-body 
combinations . 


The Munk -Jones Slender -Body Theory 


It might be well to point out here the essential differences between 
linearized theory as used in this paper and the "slender-body" theory 
initiated by Munk and R. T. Jones (refs. 1 and 2). In the M unk -Jones 
theory, which has been applied and extended by numerous authors (see, in 
particular, refs. 3 to 6 ), the first term of the Prandtl-Glauert equa- 
tion (12) is neglected in comparison with the second and third. The rate 
of change of the streamwise perturbation velocity is considered to be 
small in comparison with the rate of change of the vertical and lateral 
perturbation velocities. The application of slender-body theory involves 
solutions which satisfy Laplace's equation in planes normal to the free- 
stream direction, and the streamwise coordinate enters only as a parameter 
introduced by the body or wing-body cross-sectional area. For an a.x~t al 
distribution of singularities the slender-body counterpart of equation ( 19 ) 
is: 


4rtcp(x,r) = 2S b '(x)ln — - Sfc' (x)Zn[x(z-x) ] - f ^ ( x ) ( x ) ^ 

2 4 , I x-x| 


a^xjcos n 9 b n (x)sin n 0 


Vj ^rLWcc 

n=i 


sin n 6 ~ 
n 
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where St,' is the first derivative of the body cross-sectional area and. 
the parameters a n (x) and b n (x) represent the multipole strengths. 

The above relationship is a valid approximation only in the near 
vicinity of the body. It should be noted that the slender-body potential 
for a multipole of order n = 1 or greater depends only on the parameters 
a n (x) or b n (x) and does not involve a line integration as in linearized 
theory. The Mach number effect enters only through the first term. 


THE NONLIFTING WING IN SUBSONIC FLOW 


The nonlifting wing in subsonic flow may be represented by a suitable 
planar distribution of three-dimensional sources according to equa- 
tion (13). If the wing leading and trailing edges are straight-line 
elements, and if the wing profile is geometrically similar at all spanvise 
stations, the evaluation of the flow over the wing can be best performed 
by a consideration of elementary source filaments. In this way, the 
double integration indicated by equation ( 13 ) is reduced to a single 
integration. 


Untapered Wings 


The flow about thin untapered wings has been studied in considerable 
detail (see, e.g., refs. 10 to 13 ). The following discussion presents 
certain additional considerations useful for numerical techniques. 

The elementary source filament in incompressible flow .- Consider the 
source filament of constant strength E/Uoo per unit length (incompressible 
flow) shown in sketch (e). The perturbation velocity Sv (all velocity 





Sketch (f) 








18 


NACA TN 3964 


normal and tangential velocity components induced at the point (x,y,z) 
are 


v (n)( X,y ' z ) 

V( t j(x,y,z) 

The velocity components in the 

V (x) = V (t) 

v (y) = v (t) 

V ( Z ) * %) 


= E 

= 4rtl^N 

_ E 

Wtyr 

x,y,z 
sin A-v 
cos A+v 
sin w 


(sin t^+sin p 2 ) 


' (cos |i 2 -cos 1 ^) 

directions are 

. .cos a) cos A 
(n) 

, .cos w sin A 


(38) 

(39) 


(ho) 

(hi) 

m 


The sweptback wing . - The relationships given by equations ( 38 ) to (42 ) 
can be used to determine the compressible flow over thin, nonlifting, unta- 
pered wings if the Gothert rule is applied and if the source strength is 
defined as 


-i- = 2 jcos Ag^Z’dx (43) 


where Z* is the x-wise slope of the airfoil surface. Since only 
untapered wings are being considered here, the variable x is to be 
replaced by x-x and an integration performed from x = 0 to x = c. 

The factor cos Ag appears in equation (43) since it is more convenient 
to perform the Integration in the streamwise direction rather than in the 
direction normal to the oblique source filament. The absolute value signs 
are used in order that the source strength shall be properly defined for 
left-hand wing panels as well as for right-hand wing panels. 
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Consider idle sweptback wing skown in sketch, (g) where a coordinate 
system made dimensionless with respect to the root chord is used. The 


Uco 



axial, lateral, and vertical perturbation velocities for a compressible 
flow are to be obtained by evaluation of the following: 


Right-hand wing Left-hand wing 


2*3 2 u(l,il ,5) = |cosA e[jT 3Z«f (x) d|+jcos Ae| J 3Z'f (x) d| (44) 

o 0 

2rtPv(S,ri,£) = |cos Ae| J PZ'f , jd|+|coB | J P z ' f ( y ) d i" (^5) 

o o 

2*0v(S,ti,O = [cos Ae|J" pz'f ( z) dT+| cos J PZ'f^jd | (46) 


The first integrals on the right-hand side of equations (44) to (46) 
represent the source fil aments on the right-hand wing panel and the second 
integrals those of the left-hand wing panel. The quantities £ f x \> £ ( Y )’ 
f( z ), for the right-hand wing are: 

kcos u -cos /sin p +sin p n 

f (x) = ^ i ) sin n J cos " cos ^ 
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cos 


u -cos 

_2 

N 


n. 


^ cos Ae+(^- 


sin jj. +sin 


N 


cos co sin Ag 


(48) 


r (z) 


sin u. +sin p 
*1 ^2 


H 


l ) 


sin co 


(49) 


Although these quantities can be expressed in algebraic terms with 
the help of sketch (f) and then substituted into the expressions for the 
velocity components (eqs. (44) to (46)), it is more convenient for numer- 
ical techniques to retain the identity of certain elementary quantities 
because of the extensive repetition involved in obtaining the final 
velocity components u, v, and w. For the ri^it-hand panel of a sweptback 
wing, the following are required: 


N = 7(fH) 2+R s 2 


co = tan 


-1 K 
*2 


-1 

Li = ten — 

N 


ll = tan’ 
2 


1M 






- R, 


cos Ag 1 


H 


(50) 


where 


R 1 = Pn cos Ae+(|-|)sin Ag 

(51) 

R e = Ptj Bin Ag-(|-F)cos Ag 

(52) 


(It should be noted that in the plane of the wing N reduces to N = |R 2 j 
and cos co = Rg/| | = sign B z . ) 
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The relationships given 
"by equations (50) to (52) were 
derived for a right-hand, 
sweptback wing panel. These 
results apply to other wing 
panels if the proper sign is 
affixed to cos Ag and 
sin Ag according to sketch 
(h). The distance p(b/ 2 ) 
is negative for left-hand 
wing panels. 

The infinite oblique 
wing . - Formulas for the per- 
turbation velocities induced 
by the infinite oblique wing 
can be obtained by use of the 
previous relationships if the 
infinite oblique wing is considered 
to be composed of two semi -infinite 
wing panels as shown in sketch (i ) . 

For the left-hand wing panel, 

= it/ 2 , R 1= :-Pt] cos A e -(|-I)sin Ag, 
and Rg = -Pq sin A e +(|-f)cos A e . For 
the right-hand wing panel, p £ = jt/2, 
R 1 = Pq cos Ag+(£;-F)sin Ag, and 
Eg = Pq sin Ag-(|-|)cos Ag. The use 
of these definitions together with 
the definitions for R, p , and u 
given in equation (50) lead to the 
following equations for perturbation 
velocities : 


Ueo 



$TJ 


Sketch (h) 



■ Pi) 


Left-hand 
wing panel 


Right-hand 
wing panel 


Sketch (i) 


P 2u (i,n,£) = - 


cos 2 Ag J 1 [Pq sin A e -(|-T)cos A e ]pZ'd| 


o (P?) 2 +[Pt] sin A -(g-g)cos A _] 2 


Pv(£,q,£) = P tan AgU(|,Ti,0 


cos AgP£ p 

= - J 


pZ'd| 


o (P£) +LPt] sin Ag-U-^cos A e ]‘ 
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Since the origin for q is immaterial for the infinite wing it is more 
convenient to write these equations for = 0, 


pv(S,£) 

Pv(|,£) 


cos Ae r 1 (g-l)pz’dl 

* o (T-ft) a +0£/cofl A e ) 2 

13 tan A e u(£,£) 

n P 1 PZ'dj 

* cos A e J 0 (l-I) £ + 0^/cos A e f 


(53) 

(54) 

(55) 


where sin Ag and cos A g are defined hy equation (8). 

When the perturbation velocities are evaluated in the plane of the 
wing (£ = 0), it is permissible to set £ = 0 in equation (53) hut not in 
equation (55). However, equation (55) reduces to w = Z' when evaluated 
in the plane of the wing. 

If equation (53) is evaluated in the plane of the wing. 


u(S,0) 



Z*d| 

I-l 


(56) 


it is seen that the pressure coefficient for thin -airfoil theory (C p = -2u) 
c han ges with Mach number according to the factor 


_1_ = 1 = 1 
^ ,,/sin 2 A+P 2 cos s A ^l-M 2 cos 2 A 


(57) 


For the unswept wing this factor reduces to the Prandtl-Glauert factor 
1/p. The dependence of the compressibility effect on the Mach number 
component normal to the sweep line was pointed out by R. T. Jones, 
reference 14. 

The leading-edge singularity for round-nose airfoils . - For round-nose 
airfoils the slope Z' is infinite at the nose and changes rapidly in 
value in the imm ediate vicinity of the nose.* For numerical work it has 
been found convenient to replace the airfoil slope Z 1 by 
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z* = 


and the source strength (eq« (4-3)) by 


(58) 


jf- = 2 1 cos A e |pz' m dx 

00 


(59) 


This modification is illustrated in 
sketch (j) for an KA.CA 64A008 profile. The 
effect of this modification is to reduce some- 
what the absolute values of the calculated 
perturbation velocities in the immediate vicin- 
ity fore and aft of the nose. In reference 15 
it is suggested that the following correction 
(due to Riegel) be used. 


U(x,Z) 


Ul(*>o) 

n/i+(Z«) 2 



Sketch (j) 


where U is the total velocity over the airfoil, the subscript L refers 
to the linearized value. This is equivalent to 

U L + 1- n/i+(Z ! ) a 

7i+(z') 2 7i+(z’) 2 

which has the liim.ting value u = -1 at the nose. However, it is not 
possible to apply the Riegel rule at locations slightly forward of the 
airfoil nose. 


The use of the modified airfoil slope Z' m permits one to calculate 
reasonable and consistent perturbation velocities everywhere except at the 
wing leading and trailing edges where first-order theory gives infinite 
values except when the wing profile is cusped. For numerical work the 
calculated streamwise and lateral perturbation velocities in the plane of 
the wing have been term i nated with values -cos 2 A and cos A sin A, respec- 
tively, at the wing leading and trailing edges. The pressure coefficient 
may be estimated by substitution of the values u = -cos 2 A and 
v = cos A sin A into equation (9). 
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The inherent singularity at C - 0 as H -» 0 .- If a numerical 
technique Is used, the inherent singularity for £ = 0 at K 0 can 
he circumvented in the following manner. For the inf ini te oblique wing 
the singularity at £ = 0 as f g is avoided by writing equation ( 53 ) 
in the form 


P 2 u(^0) 


cos Aq 

it 



wag r s 0>z'> 8+ s-(PZ')g.E 
T-t + 4> S 


d E 


XSSL 


( 60 ) 



Sketch (k) 


where 5 is any convenient 
subinterval bracketing the 
singular point, see sketch (k) . 
The integrand of the second 
integral of equation ( 60 ) has 
the following finite value at 
the singular point 


lim 
&-> o 


(PZO^-O 2 ')^ 

5 


- 2pZ"(&) 


where Z"(g) is the second 
derivative of the airfoil ordi- 
nate evaluated at the singular 
point. 


For swept wings an inherent singularity in the plane of the wing 
occurs when | has the value 

T* = £-(3t] tan A^ ( 6 l) 


The singular integrals of equations (44) and (45) for £ = 0 can be avoided 
by writing these integrals in the form 


r 
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[cos A pZ’f^^dl 


*elj *“ "(x) 

o 


n£ "S n 1 __ 

-I = j cos Ag[ J pZ*f( x )d|+|c°s Agl !J + 


l*+S 


s 

jcos A e | / ■ [pz** (x) ] +[pz'f( x) ] 


_* _ 
I -5 


- 

I +8 


dS 


(62) 


_*• 

>£ -S 


| cos A e [jf pz‘f (y) d| = | cos AgjJ" pZ l f( y) d|+jcos A e \f P z ' f ( y )^ + 


r+s 


i“ A /' [ pz,f ( y )j , + [»’*«!. 

I 1-5 | +6 J 


dS 


( 63 ) 


where 6 is any convenient svbinterval bracketing the singular point 
= (where s’ = 0) . The integrands of the third integrals on the 
right-hand sides of equations ( 62 ) and ( 63 ) have the limiting values 


and 


2 


2 


sin AgPZ’ 


cos AgpZ 1 


1 1 


3(5/2) 


3t] 

cos Ag 

cos Ag 

cos Ag 


1 

1 

3(5/2) 

3t] 

3*1 

cos Ag 

cos A e 

cos A e< 


- 4-3 Z" 


+ 4 tan AgPZ" 


respectively, at the singular point. (This technique for avoiding the 
singularity in the plane of the wing was illustrated previously for the 
infinite oblique wing in sketch (k).) 
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The distortion velocity at the center Une of a sveptback -wing . - At 
the center li ne of the sveptback wing a unique result is obtained -whereby 
the streamwise velocity is found to be equal to that for the infinite 
oblique wing plus a distortion velocity of the Ackeret type (proportional 
to Z 1 ), see reference 11, 12, or 13 . This result is easily obtained in 
the following manner. First consider the right-hand, sveptback wing panel. 
The streamwise perturbation velocity at the center line can be obtained by 
use of the relationships of equations (44) to (52) with £ = 0, ^ = it/2, 

U = Jib, |. There is obtained 


2jtp 2 u( 1,0,0) 


lim 

Pt)-=*o 



pz'df 

Pn si 11 A e -U-T)cos A e 


cos 




frlPZ'dg 


0 [0T1 sin Aq-U-Dcos Ae]7(|-F) 2 +(Pil) 2 J 


(64) 


In the first integral q can be set equal to zero and the result is equal 
to one half that for the infinite oblique wing. The results for the left- 
hand, sveptback wing panel can be expressed in a form identical with 
equation (64). Thus, the distortion velocity at the wing center line, 

(Ptl = 0), 


u d u center ” u oblique . 
line wing 


is given by 


P u (l) = Hm 

pq-»o 


cos Ag J' 


PnPZ'dl 


[Pq sin A e -(g-|)cos J 


It is clear that the above integrand is zero everywhere for 3q = 0 
except at the singular point where the integrand is O/O. In order to 
obtain the limit, an integration is performed from |-Pq tan Ae-5 to 
£-Pq tan A e -€ and from g-pq tan Ag+e to g~j3q tan A e +5 with PZ’(^) 
replaced by ' (3Z 1 (|)+PZ”(|)J. The Cauchy principal value is obtained by 
letting e -> 0 and the final result is independent of 5, 
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u d(5) = - 


cos AqZ’^) 1 +sin Ag 

1 Zn ■■■ — — 

1 -sin Ag 


(65) 


It should he pointed out here that, in the original body contouring method 
of Kiichemann (ref. 13), the distortion velocity at the wing-body Juncture 
is assumed to be that given by equation ( 65 ). 


For sweptback wings of finite aspect ratio the perturbation velocity 
at the center line for £ = 0 is 


|cos 

u = u oblique + u d + 


wing 


cos A_ 




1 - 


1 l B 


st3‘ 


pcyg) R 

cos Ag i 


a p(b/2) 

7_cos Ag 



sin A„ 


Pfo/g) 

(_cos Ae 


3Z*dg 


12 


- S, 


where R x = (|-|)sin Ag and R 2 = (|-J_)cos Ag. (The above integral is not 
singular at the singular point f" = g * ) 

In the immediate vicinity of the tip of a sweptback wing the 
streamwise perturbation velocity may be approximated by the following: 


u = u , . - — 

oblique 2 

wing 


Wings Having Tapered Plan Forms 


Ihe thin, nonlifting wing with tapered plan form has been considered 
In references 16 and 17 where an evaluation of equations (l4) to (l 6 ) was 
performed by integration in the streamwise direction first, and nu me ric al 
results could be obtained only by laborious expansions into inf ini te 
series. However, as will be shown by the following analysis. It is pos- 
sible to compute the flow about the tapered wing by use of the source 
filament concept in a manner similar to that for the untapered wing. 
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Fundamental relationships . - Consider the tapered wing panel in 
sketch (2 ) (the equivalent plan form according to the affine transforma- 
tion of the GSthert rule is shown), where the wing profile is geometri- 
cally similar at all spanwise stations j that is, Z' is constant for a 
given angle a. The integration over the wing plan form indicated by 



Sketch (l) 

equations (l4) to (l6) is obtained by first integrating with respect to 
s', keeping the angle a constant. The results of this line integration 
may be thought of as representing the influence function for a line, or 
filament, of source-like singularities. 

The sources lying within the tapered, shaded area of sketch ( l ) are 
replaced by an equivalent source line, or filament, having strength 

2 Ksiew - T ) to 
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V 


per unit length of s', 
see sketch, (m) . The 
perturbation velocity 
components induced hy 
this source filament 
in directions normal 
and tangential to the 
filament are 



Sketch (m) 




and the streamwise, lateral, and vertical velocity components induced by 
the filament are given by equations (40) to (*i2). The quantities N, 

R£, u, and are the same quantities defined previously by equa- 
tions ( 50 ) to ( 52 ). The parameter see sketch (l) , is related to the 
angle a by the following relationship 
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a 


<• 


COfl cot a„ - 
LE 



(68) 


For numerical work it is convenient to use the parameter | as the 
variable of integration; hence, da in equations (66) and (67) is to he 
replaced by- 


da 



l cos A I _ 

e_i_ ^ 

P L/cos Ag 


(69) 


The sweptback wing .- The perturbation velocities induced by the 
right-hand wing panel shown in sketch (l) are obtained by evaluation 
of the following: 


2rt0 a u jcos Ae|PZ>rg (x)+ h (x) ld| 


(70) 


2jtPv 


= / l COE Ae| 


0Z 1 


_ g (y) +il (y)J dI 


(71) 


2rtPw 


= f ( cos Ag| 


PZ’ 


S (z) +ll ( Z )J dI 


(72) 


where 


6 (x)=^- 


g (y) 


s 


■0 


(zrv 1 


*1 

_\ 

"/cos 

M-g “COS h x n 

pL /cos 

V 

y 

N ./ 

R i 

_\ 

■/cos 

}i 2 -cos 

p L/cos 

A*/ 

y — 

N , 

A 

\ 

/sin 

u +sin p \ 
2 1 

P L/cos 

V 

V"~ 

N / 


-^) sin Ag - (■ 


sin p 2 +sin p. - 


N 


y cos w cos Ag 


(73) 


''sin a +sin l\ 

~ — J cos a) sin Ag 


N 


(74) 


sin w 


(75) 
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l (x) 


cos a. L LV C0S / V =08 u 2 J\ 


P L /cos Ag 

(sin p 2 +sin p^sin Ag-(cos Hg-cos ji-^cos t») cos Agj- ( 76 ) 


L (y) 


i Loe Ag -A a£\] . 

M . I ^ LV =° s Hj. / v = OB i* 2 /J 


P L/cos Ag 


(sin p 2 +sin p. 1 )cos Ag+(cos p £ -cos p^cos u sin A e "j- 


(77) 


h (z) - 


P L/cos A e 


j^(cos p 2 -cos p^sin u 


(78) 


When the parameter pL/cos Ag becomes inf i n i tely large the h 
functions vanish and the g functions reduce to the relationships given 
previously for the untapered wing (eqs. (47) to (49)). If the plan-form 
convergence factor, defined as 


K a 


P 


C R~ C T 

P(b/2) 


(79) 


is small, say less than 0.5, the h functions may be neglected in 
equations ( 70 ) to ( 72 ). 

Evaluation of the streamwise and lateral perturbation velocities in 
the plane of the wing .- If a numerical techn i que is used, the singularity 
in the plane of the wing at |* (that is, at N = 0) can be avoided by a 
procedure similar to that used for the un tapered wing, with the exception 
that the log term appearing in h£ x ^ and ^(y) mus "^ evaluated sepa- 
rately. Eor £ = 0 it is convenient to write equations ( 70 ) and ( 71 ) in 
the form 


32 


MCA TN 3964 



■where 8 is any small, convenient subinterval (say 5 = 0 . 05 ) "bracketing 
the singular point. 

The first three integrals of equations (80) and (8l) are not singular. 
The third integrals of equations (80) and (8l) have the following limiting 
values, as 8 0 for a right-hand wing panel. 
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sin Agcos Ae^Z 1 ^ pL 


PL/cos Ae VC ° S ^ 


_J±\ 

cos A q ) 


1 

P(Tj/2) _ g4 

cos Aq cos Aq 


1 

3g 

cos Ag 


and 


4 


P L /cos Ag 
cos Ag 



cos 2 A e PZ 1 / pL 
pL/cos Ag Vcos Ag 



1 

Pfo/2) _ 34 

cos Ag cos Ag 


Pq 

cos Ag 


4 


tan Ag(P L/cos Ag) 
cos A e 


cos e A e PZ’ 
P L/cos 


df 


respectively, at the singular point. (The indicated partial derivatives 
may he evaluated graphically. ) 

The integrands of the fourth_integrals of equations (80) and (8l) 
are logarithmically singular at £ — > |* and do not change sign at the 
si ngula r point I*. These integrals can he evaluated approximately if 

P(h/2) Pq 
cos Ag ~ cos Ag 

sin A e 


Pq/cos A^ 

< , 


sin Ag 


as follows : 
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I 


£*-& 


|*+5 


cos AgPZ’h^dl 


4& sin AgCOs AgPZ* 
pL /cos Ag 


In < 


/ Pn 

[P(V2) PH 

/ cos a®. 

_cos A e cos Ag_ 


C 1 -^-) 5003 ^ 


I 


£*-8 


r*+s 


cos AgPZ'li/ ,4| 


45 cos 2 A e PZ' 
P L/cos Ag 


In < 


/ PtJ 

p(b/2) 

Ptl 1 

cos Ae 

_cos Ag 

cos Ag_ 


1 


( 1 -^) 6 003 A e 


At the wing center line the lateral velocity is zero "because of 
symmetry. For numerical work, the streamwise perturbation velocity can- 
not be obtained from equation (80) by setting Ptj = 0 since the function 
g^ contains an indeterminant form for Pq = 0. However, equation (80) 

can be used at the wing center line if the quantity 2rtP 2 u^ is added to 
the right-hand side. At the wing center line the third integral of equa- 
tion (80) has the following limiting value as 5 -3>- 0 : 


4 


cos Ag 


sin A e PZ' 
,p(b/2)/cos A e 



At the wing center line the fourth integral of equation (80) contributes 
the approximate value, 
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46 sin AgCos AgPZ’ 
P L/cos Aq 


In 


0O>/2) 

2 7— 

cos Ag 
6 cos Ag 


The above results for the right-hand, sweptback 'wing panel can be 
used for other wing panels if the signs of cos Ag and sin A@ are changed 
in accordance with the sign convention of sketch (h). The angle cr, which 
is the complement of the Gothert equivalent sweep angle Ag, is considered 
to be a positive angle (in the second quadrant) when measured counterclock- 
wise as in sketch (z). For the left-hand sweptback wing panel, the angle 
cr is measured clockwise and hence is a negative angle (in the fourth 
quadrant). The distances p(b/2) and pL are considered to be negative 
for left-hand wing panels. 


THE SUBSONIC FLOW ABOUT NONLIFTING BODIES 


The linearized spatial flow about bodies can be represented by 
suitable axial distributions of sources and multipoles according to the 
general solution ( 19 ) where the body axis and the axes of the multipoles 
are aligned with the wind axis. It Is assumed that the body cross- 
sectional area distribution terminates with zero first derivative, 

S'(0) = S*(z) = 0 . It is further assumed that the function describing 
the body shape, R(x,0), Is single-valued in 0 and possesses smooth and 
continuous partial derivatives. 

When equations ( 32 ) to (34) are used in conjunction with the boundary 
condition (6) there results the following general requirement to be sat- 
isfied in the evaluation of the strength distribution functions A a and 

BjlI 



K n 

4*(PR) I1+1 


[A n (x)sin n0-B n (x)cos n0] 


5 PR 
PRd0 


I 


n=o 


Sn 

4jt(PR) n+i 


[Aj^xJcos n0+B n (x)sin n0] 


= 0 


( 82 ) 
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Axially Symmetric Bodies 


Bor circular todies , the source solution (n = 0) is the only solution 
involved and the boundary condition (82) reduces to 


A q (x) = 2*3R &£■ = P s S^(x) 
dx 


(83) 


■where S'-^x) is the first derivative of the body cross-sectional area. 
The perturbation-velocity components are 


u(x ’ r) ' t I 


(x-x)S' b (x) 

4ot ^ r ( —\2 q 2 2-, 3 ^ 2 
o [(x-x) +p r ] 


dx 


(84) 


T (r) 


M ’ £ I 


SV(x) 


'O [(x-5f 


dx 


(85) 


and the pressure coefficient is to be calculated from 



-2u-v 


(r) 


2 


When the strength distribution of an axial source line cannot be 
expressed analytically, the integration of equations (84) and ( 85 ) must 
be performed by a suitable numerical means. Although the integrand is 
not singular at x = 0 for either equation, the integrand of equation (84) 
changes sign and varies greatly in magnitude near x = x for small Pr. 
Consequently, the numerical evaluation of equation (84) is simplified 
considerably if the region of integration is separated into subintervals 
(see sketch (n)) and this equation written in the form 

x < 1/2 


(0,0) (X,0> (2x,0) 


( 1 , 0 ) 


x > 1/2 


( 0 , 0 ) 


(2x-l,0) (x,0) 

Sketch (n) 


( 1 , 0 ) 
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p2X — r 

U(X ’ r) ' t 1 




O E(\r+Pr ] 


_1_ 

4it 


/ 


(x-x)S’ b (x) 

w r , _,2 _2 2.3/2 
Z-2X L (x-x) +p r ] 


dx 


> (86) 


x<i 

2 


u(x,r) 


p2—x 

-if - 

4jt vJ , 


1l jr J r/ — .2 2 2 .3/ 2 

^ o C(S X ) +P r ] 


S’ . (x-Bj-Sy**^) 


d^ + 


i r 

w J 


:x " 1 (ac-XjS^C*) 

r / — \2 -2 2 .3/2 

o [(x-x) +p r J 


dx 


> (87) 


x > |- 


■where 6^^ = x-x. 


Nonsymmetrlc Bodies 


For the flow about bodies which do not possess axial symmetry, a 
consideration of both sources and multipoles is required. The body 
shape may be designated by 
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R(x,0) = R(x) + 



n=i 


( 88 ) 


■where R is the radius of a circular tody having an identical axial 
distribution of cross-sectional area S^x) and -where £R n is the incre- 
mental change in tody radius associated with the nth order multipole 
distribution. 

In most practical cases the tody will possess vertical and lateral 
symmetry and will te nearly circular in ctosb section. In this case the 
boundary condition (82) can te simplified to the following. 


n 



jrrr A n (x)cos n0 = 0 

(SR) + 


(89) 


In order to solve equation (89) for the strength functions An(x) it 
is convenient to express the tody slope in the form (primes denote 
differentiation with respect to x). 


n 

V -1 


PR‘(x,0) = pfi’(x) +2^ 


PAR , n (x,0)(PR) 


n+i 


(PH) 


n+i 


(90) 


It is assumed that the even function PAR' = p(R , -R* ) can te 
represented by a Fourier cosine series and thus equation (90) becomes 


PR 


’ (x,0) = pfi*(x) 

n=2,4,e , . . . 


a n (x)cos n0 

(pR) n+1 


(91) 


where 


k r* /£ 

a n (x) / PAR'(pR) n+1 cos n0 d0 


( 92 ) 
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If the relationships (91) and (92) are substituted into ( 89 ) and if 
the coefficients of cos n0 are set equal to zero, there results in the 
following solutions for the source and multipole strengths: 


Aq = 


1 


An 


= ±i r 

Ku 1 


PAR 1 ((3R) n+1 cos n0 d0 


(93) 


where 


AR’(x,0) = R'(x,6)-B'(x) 


Once the source and multipole strengths have been determined, the 
perturbation velocities may be evaluated by the use of equations (20) 
to ( 31 ). 


WING-BODY COMBINATIONS 


It is assumed that the 
wing-body combination pos- 
sesses both lateral and ver- 
tical symmetry and a coordi- 
nate system is used which is 
dimensionless with respect 
to the junction chord and 
with origin as shown in 
sketch (o). The body radius 
is denoted by R and S is 
used, when the body is non- 
circular, to denote the 
radius of an equivalent cir- 
cular body having identical 
cross-sectional area. 

The boundary condition 
to be satisfied at the body 
surface is 
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&R SR . SR 

S5j R 


sr + _sr_ y 
s© RSe Li 


v (e) -V (r) " V (r) 
n v 


n 

'X v (r) -° (54) 

n=o Aj v 'n 


where 


and v, 


are the radial and azimuthal velocities induced 


(r) v “ v (e)^ 

by the wing (which is assumed to extend through the body from wing tip 
to wing tip as indicated in sketch (o)). The radial velocity induced by 


the axial sources is denoted by v/ \ , 

n=o 

the contributions of the axial multipoles. 


and 


r (r) 


and v 


n 


(9) 


represent 


n 


The strengths of the axial sources are chosen so as to provide the 
required body volume or, more specifically, the body cross-sectional areas, 
and to e limin ate the average of the wing-induced radial velocities at each 
axial station. The average wing-induced radial velocity is defined as 


r 

2 / 

\ y M 

* J 

L W- 

av 0 


,J c/2 


m 


( 95 ) 


The strengths of the axial multipoles are chosen so as to complete 
the cancellation of the wing-induced radial velocities and, for non- 
symmetric bodies, to provide, in addition, the required distortion 
AR(|j,0) = R(£j,0)-R(£j ) from the circular shape having equal area. 


Determination of Axial Source and 
Multipole Strengths 


Bodies having circular cross sections . - When the body shape is 
circular in cross section, the boundary condition at the body surface is 


0 


3R 



-I 

n =n 



( 96 ) 


where the radial velocities, v 


Wn’ 


induced by the axial sources and 


multipoles can be expressed in terms of the strength functions A n (^.) 
by the use of equations (32) and (3I*). “ 
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Since we are considering only symmetrical, no nli fting wings, it is 
possible to represent the wing-induced, radial velocity by a Fourier cosine 
series involving only the even harmonics, 

Pv (r) = 20 +a 4 (^)cos k8 + . . . ( 97 ) 


where 



(98) 

(99) 


thus. 




(r) 


= P 


w 


( r ) 


1 k T' r* /z 

+ — } cos n 0 J PA v^ r ) cos n 0 d 0 ( 100 ) 

wJ av ^, 4 , 6 ,,.. ° w 


A consideration of equations ( 32 ), ( 96 ), and (100) leads to the 
following definitions for the axial source and nrultipole strengths : 1 


A q = 2*!® 


{ p £ ■ M J 


( 101 ) 


An 

m n 



(102) 


Bodies having noncircular cross sections .- If the body shape is 
noncircular (but possesses both lateral and vertical symmetry), it is 
convenient to represent the body shape by 

R(Sj,e) = R(ij)-^a( 6 j, 0 ) ( 103 ) 

Anhe quantity Av/ % appearing in equations (100) and (102) 
v 

replaced by v/ \ if so desired. 


can be 
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■where AR is an even function in 0. The axial source strengths are 
given by 


\ T( r ) 

o 



{s'b - 2/^ 



(104) 


~2 

where = jtR is the equivalent body having equal cross-sectional area. 
The strengths of the axial multipoles are to be chosen so as to complete 
the cancellation of the wing-induced radial velocities at the body surface 
and to provide the required distortion AR from the equivalent circular 
shape. 


The boundary condition 


for determining the multipole strengths is 


0 


(Bv/ \ 

5AR SR 

atj R S0 





n=a,-t 


>••• 


03B) n+1 


(' 


cos n0 


— sin n0 ) = 0 ( 105 ) 

R&0 / 


If the above equation is evaluated for a constant f y at appropriate 

azimuthal angles 0, a set of simultaneous equations is obtained involving 
the unknown parameters If the body shape is nearly circular the 

azimuthal derivative <3R/c 50 may be neglected, and, if the body slope is 
approximated according to equation (90) , the multipole strengths are given 
approximately by 


•^n _ 16 

(PR) n = 


r 


|3AR(pR) n+1 coB n0 d0 - 


16PR 

K n 




cos n0 d0 
w 


(106) 
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Pressure Distributions 


Once the source and nrultipole strength distribution functions have 
been dete rmin ed the pressure distributions may be calculated by first 
obtaining the perturbation velocities and then using equation (10) or (ll). 
At the body surface the pressure coefficient is approximately 


°p - + i ”-•) - (0 

n=o J 

while at the wing surface the thin-airfoil approximation can be used. 



n=o 


The general solutions for the perturbation velocity potential 
(eqs. ( 13 ) and ( 19 )) behave properly at large distances from the wing 
or body, that is. 


u,v,w 0 

as the distance becomes infinitely large. The theory presented in this 
report can be used to calculate pressure distributions in the neighboring 
flow field as well as at the wing-body surface. 


ILLUSTRATIVE EXAMPLES AND COMPARISON WITH EXPERIMENT 


In order to illustrate and, to a certain extent, justify the practical 
use of the present method for calculating pressure distributions for wing- 
body combinations, a comparison of some theoretical and experimental 
results will be presented. Because of the current interest in area-rule 
concepts, a comparison of theoretical and measured wing-body juncture 
pressures will be made for two sweptback wings of aspect ratios 3 811 d. 6 
in combination with basic Sears-Haack bodies of revolution and in combi- 
nation with the basic bodies indented according to Whitcomb’s transonic 
area rule (ref. 18). The experimental data used in the following com- 
parison with theory were obtained in the Ames 14-foot transonic wind 
tunnel. 
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Theoretical Pressure Distributions at = O .85 for 
an. Aspect-Ratio-6 Wing in Combination With a 
Basic and Area-Rule Indented Body 


The plan-form details of the aspect-ratio-6 wing in combination with 
a basic Sears -Haack body of fineness ratio 11 and in combination with the 
basic body indented according to the transonic area rule are presented in 
figure 2. The procedure for calculating the pressure distributions at 
= 0.85 will be described in detail for the indented body configuration. 

Perturbation velocity components induced by the wing alone . - For this 
wing of large aspect ratio the taper effect is small, and the formulas 
given previously for determining the flow over untapered wings may be used. 
The flow over the wing at = O .85 was obtained by the numerical evalua- 
tion of equations (44) to (46) using the relationships of equations (47) 
to (52). (A constant sweepback angle of A = 40° was used. ) 

The locations of isobars at the wing surface are shown in figure 3* 
These results were obtained by use of the thin-airfoil approximation 
Cp = -2u. The calculated perturbation velocities at the indented body 
location are presented in figure 4 for several azimuthal locations (see 
sketch (o)). The radial component of the perturbation velocity at the 
body surface, 



= v cos 0 + w sin 0 


(107) 


and the average radial velocities, (eq. (95) ), are presented in figure 5* 

Determination of the axial source and multipole strength 
distributions . - For this configuration the body shape Is circular and 
the axial source and multipole strengths were obtained by the use of 
equations (101) and (102). The variation of source strength with axial 
station is shown in figure 6. Equation (102) was evaluated numerically 
(see fig. 7 for typical variations of 3A v ( r ) cos n0 with 0) and the 

results for multipoles of order n = 2, 4, and 6 are presented in 
figure 8. 

Wing-body juncture pressure distributions .- The streamwise 
perturbation velocities induced by the wing planar sources, axial sources, 
and axial multipoles of order n =2, 4, and 6 _for the aspect-ratio-6 wing 
in combination with the indented and basic bodies are shown in figures 9 
and 10. The difference in the flow fields about the two configurations is 
primarily the result of differences in the strength distribution of axial 
sources (see also fig. 6). The Influence of the axial multipoles on the 
flow field Is comparatively small and decreases rapidly with an increase 
In the order n of the multipole . 


F 
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A direct comparison of theoretical -wing-body juncture pressure 
distributions for the basic and indented bodies is presented in figure 11 . 
Also included in figure 11 is the theoretical pressure distribution at the 
wing ml d-s emi span (which is sufficiently far removed from the model center 
line and wing tip so as to be essentially free of end effects in subsonic 
flows-). For the indented body the juncture pressure distribution is almost 
identical with the pressure distribution of the mid-semispan. Additional 
calculations at other spamwise stations indicate that the isobar pattern 
on the wing for the indented combination (excluding the local regions near 
the wing tips) Is essentia lly two-dimensional, that Is, the Isobars tend 
to follow the local sweep lin es of the wing. 


Comparison With Experiment 


A comparison of some theoretical and measured pressure coefficients 
at Moo = 0.85 for the aspect-ratio -6 wing in combination with the indented 
and basic bodies is presented in figures 12 and 13 . The comparison of 
theory and experiment is made for two body azimuthal locations , 2 0 = 0 ° 
and 0 = 90 °. Hie measured pressures at a row of wing orifices In the 
immediate vicinity of the body (see fig. 2) are also Included In fig- 
ures 12 and 13 since it was felt that a comparison of theory and experi- 
ment at the wing-body junction would suffer somewhat from boundary-layer 
effects in the corners of the wing-body juncture. In general, the theory 
and experiment are in reasonably good agreement. 

Similar comparisons of theory and experiment for the aspect-ratio-3 
wing and body combinations shown In figure 14 are presented in figures 15 
and 16. For this case the perturbation velocities induced by the wing 
alone were evaluated by numerical means using the theory for wings of 
tapered plan form. The location of isobars at the wing surface are shown 
in figure 17. A convenient formula for obtaining the variation of the 
local sweepback angle with junction-chord station % is 

«J 


tan A e (ij) = tan A e (0) (108) 


2 At the wing-body juncture the pressure orifices along the sides of 
the bodies were located as close as possible to the upper surface of the 
wing. The designation 6=0° for the side row of orifices is used here 
for convenience, although, at the wing-body juncture, the or ifi ce loca- 
tions differed slightly from 0 = 0 ° because of the finite thickness of 
the wing. 
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CONCLUDING REMARKS 


In this report a theoretical method haB "been developed for determining 
the linearized subsonic flow about symmetrical, nonlifting wing-body com- 
binations. Particular emphasis has been placed on developing practical 
methods for use in obtaining numerical results. As part of the present 
analysis a theory was developed for determining the flow about noncircular 
bodies by the use of axi ally distributed sources and nrultipoles. The flow 
over thin, nonlifting airfoils having symmetrical profiles was also studied 
in considerable detail and it was shown that the concept of the oblique 
source filament could be used for wings having tapered as well aB untapered 
plan forms. This is of considerable importance since numerical results can 
be obtained for the tapered wing almost as easily as for the untapered 
wing. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Feb. 5, 1957 
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Figure 1.- Continued. 
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Figure 5. - Wing-induced radial velocities at the area-rule body 

Mqq = 0.85. 
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Figure 6.- The variation of axial source strength with junction -chord station £ j 
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Figure 7- - T3ie variations of PAv, . and PAv/ v cos n0 with 6 







58 


NACA TN 3964 



-.2 0 .2 .4 .6 .8 i.O 1.2 

Junction-chord station, 

Figure 8.- The variations of the axial multipole strength parameters 
An/CfiR) 11 with junction-chord station 
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Figure 10.- The streamwise perturbation velocities induced at the wing- 
body juncture by the wing sources, axial sources, and axial multipoles 
for the aspect -ratio -6 wing in combination with the basic body; 

Moo = 0.85- 
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Figure 11. - A comparison of calculated pressure coefficients at the 
wing-body junctures of the aspect-ratio -6 wing in combination 
with the basic and indented bodies; = O. 85 . 
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Figure 12.- A .comparison of .theoretical . and measured pressure coefficients 
at * 0-85 for the aspect-ratio-6 wing in combination with an area 
rule indented body. 
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Figure 13.- A comparison of theoretical and measured pressure coefficients 
at 1^ = O .85 for the aspect-ratio-6 wing in combination with a Sears - 
Haack body. 
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Figure 15 • - A comparison of theoretical and measured pressure coefficients 
at Moo = 0.85 f or the aspect -ratio -3 wing in combination with an area 
rule indented body. 
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Figure 16. - A comparison of theoretical and measured pressure coefficients 
at = O .85 for the aspect-ratio-3 -wing In combination with a Sears- 
Haack. body. 
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